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Summary 

! 

A complete  theory  Is  presented  for  the  statistical  effects  of  atmos- 
pheric turbulence  on  coherent  radiation  reflected  from  a diffuse  target. 
This  study,  which  is  motivated  by  the  need  to  understand  speckle  and 
scintillation  effects  on  the  operation  of  coherent  adaptive  optical  trans- 
mitter (COAT)  systems,  constitutes  a significant  advance  in  the  field  of 
turbulence  scattering  phenomena. 

Both  phase  and  amplitude  perturbations  are  taken  into  account,  and 
the  analysis  includes  multiple  scattering  (saturation)  conditions  and 
finite  receiver  apertures.  The  development  is  free  from  certain  re- 
strictive assumptions  employed  in  previous  work,  and  yields  results  for 
the  variance  and  covariance  of  irradiance  which  lead  to  clear  physical 
insight.  It  is  found  that  the  covariance  comprises  two  additive  terms 
which  represent  respectively:  (1)  the  incoherent  scattering  mechanism 
which  is  independent  of  source  spectral  width,  and  (2)  the  coherent 
mechanism  related  to  "speckles."  A technique  is  also  developed  for 
quantitative  numerical  predictions  using  reasonable  computer  time.  Re- 
sults from  an  ongoing  program  of  experimental  verification  are  given 
which  support  the  theory. 

The  effects  of  nonideal  lasers  and  diffuse  targets  are  discussed 
briefly,  including  glints,  finite  spectral  widths,  and  truncated  target 
regions.  It  is  found  that  a small  diffuse  target  (or  region  on  a complex 
target)  can  give  rise  under  strong  turbulence  conditions  to  scintillations 
with  a large  normalized  irradiance  (3),  thus  constituting  a very  strong 
fading  mechanism. 
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I.  Introduction 

This  report  extends  the  work  described  in  Ref.  1 involving  the 
statistics  of  a scintillating  field  from  a laser- illuminated , diffuse 
target  in  turbulence.  In  particular,  the  significantly  restrictive 
assumptions  of  the  earlier  treatment  have  been  dropped,  and  clear  physical 
interpretations  of  the  resulting  complete  theory  are  presented. 

The  basic  analysis,  interpretation,  and  ancillary  theoretical  con- 
siderations are  given  in  Sec.  II,  followed  in  Sec.  Ill  by  a numerical 
technique  for  quantitative  predictions  using  moderate  computer  time. 

In  Sec.  IV,  supporting  experimental  efforts  are  described. 


1.  J.  R.  Kerr,  et  al.  Propagation  of  Multiwavelength  Laser  Radiation 

Through  Atmospheric  Turbulence,  RADC-TR-76-11 1 , April,  1976,  (A024863). 


II.  Analytical  Description 


In  See.  II-A  we  present  the  complete  analysis  of  the  basic  statistical 
properties  of  irradiance  after  reflection  from  a diffuse  target  in  at- 
mospheric turbulence.  In  II-B  we  briefly  treat  further  topics  such  as 
glints  and  other  considerations  relating  to  the  characteristics  of  real 
lasers  and  targets. 

As  in  the  preceding  report,^  the  treatment  utilizes  the  extended 
Huygens-Fresnel  principle,  and  includes  the  effects  of  tirbulence  on  both 
the  path  from  the  transmitter  to  the  target  and  that  from  the  target  to 
the  receiver  (Figure  1).  Unless  otherwise  stated,  a TEMQ0  laser  is 
assumed,  with  a perfectly  diffuse  target.  Any  assumptions  are  clearly 
stated  as  they  are  introduced;  basically  however  the  treatment  is  quite 
general  and  not  significantly  limited  by  assumptions  such  as  were  utilized 
in  previously  reported  work. 


II-A. Basic  Statistics  of  Speckles  Through  Turbulence 


1.  Mean  Irradiance 

Referring  to  Figure  1,  we  assume  that  the  source  and  target  are 
much  smaller  than  the  path  length  (L) , and  that  the  distance  between  the 
receiver  and  source  is  greater  than  the  source  size  such  that  outgoing  and 
incoming  radiation  experiences  independent  turbulent  regions.  The  latter 
limitation  is  thought  to  be  inessential  owing  to  the  diffuse  target 
characteristics,  and  is  under  further  study. 

We  write  the  source  amplitude  distribution  as 


Uo(r) 


1)  exp 
o r 


2a  2 
o 


ikr' 

2F 


where  a and  F are  the  characteristic  beam  radius 

o 

respectively.  The  field  at  the  target  is  written 


and  focal  length 
from  the  extended 


(1) 


I 

! 


In  particular,  this  applies  to  the  special  cases  of  a focused  (L  = F) 
or  collimated  (F  ■+  °»)  beam  respectively. 

The  field  at  the  receiver  is  written  by  reapplying  the  Huygens- 
Fresnel  principle  to  the  field  at  the  target: 

U'(p)  exp  (p2  - 2p-p)  + i;'2(p,p)J  dp 

(4) 

where  U' (p)  is  the  field  solution  after  reflection  from  the  target,  and 
'>2  represents  the  turbulence  effect  from  the  target  to  the  receiver . The 
mean  intensity  at  the  receiver  is  then 

2 

<I(p)>  •=  <|U(F)|2>  = jf  dpjdp2  <U'(p,)U'*(p2)> 

• exp  [|£  ((pf  - p|)  - 2p.  (P]-P2))] 

• < exp  [VF.P1)  + V''2*(f.P2)]  > (5) 

2.  R,  F.  Lutomirski  and  H.  T.  Yura,  "Propagation  of  a Finite  Optical 
Beam  in  an  Inhomogeneous  Medium,"  Applied  Optics,  _10,  1652,  July  1971. 

3.  H.  Yura,  "Mutual  Coherence  Function  of  a Finite  Cross  Section  Optical 
Beam  Propagating  in  a Turbulent  Medium,"  Applied  Optics,  11_,  1399, 

June  1972. 
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Through  the  assumption  of  a diffuse  target,  the  reflected  beam 
suffers  a random  phase  delay  from  point-to-point  over  the  target,  so  that 


<U'  (Pi)U'*(Pj)>  - <I(p,)>  «J,-p2)  . £ (6 

The  factor  W k2  requires  comment.  The  dimensions  of  the  ^-function  are 
(area)-1,  and  this  factor  reflects  the  correlation  which  actually  exists 
over  a finite  area  of  order  k-2.  The  numerical  value  4s  Is  discussed  In 
Ref.  4,  and  verified  below  by  an  Independent  approach.  Using  this  In  Eq. 
(5),  the  mean  Intensity  becomes 


<I(p)> 


■ — j dP|  «|U<;>,)|7>  < e.\p  £t:(p,P])  + l )J 


where  the  mean  exponential  term  is  unity  from  considerations  of  energy 
3 

conservation.  The  resultant  mean  intensity  at  the  receiver  is  then 
simply 


<I<P)> 


-S-  fdp 

»L2  J 


<|U(p)|?> 


To  complete  the  solution,  we  use  Eq.  (3)  with  Eq.  (8).  We  note 
that  the  structure  function  gives  us  (r  » irj-r2|) 


[*l(P. 


c*P  l«i <P,rj)  + *j*(p,r2)  > - e 


r">r  the  focused  beam,  we  then  have 


;|U(p)|2>  = U0  2 JJ  drldr2  T-  — j 

-ft»1 


ik  — — — 
l'  P ' <r  l~ r2 ) 


4.  M.J.  Bc-ran  and  C.B.  l'arrent,  Jr.,  Theory  of  Partial  Coherence,  Prentice- 
Hall,  Englcwoou  Cliffs,  N.J.,  1964. 


This  describes  a target-spot  of  size  » 


L 


ka 


or 


, , for  a , o smaller 

"o  kpo  0 ° 

respectively.  Carrying  out  the  integration  indicated  in  Eq.  (8),  involving 

the  Four  ier-Bessel  integral,  we  have  finally 

2 


<Kp)>  - U 


a 

2 o 


Vj 

o 

U 2 a 2 
0 0 


-r2/a 


(12) 


o iy  <“> 

The  result  for  the  collimated  beam  is  identical,  and  in  fact 
could  be  deduced  lor  an  arbitrary  beam  focus  through  conservation  of 
energy,  as  follows.  For  Lambertian  scattering  into  a hemisphere,  the 
on-axis  irradiance  at  a distance  L from  the  target  is 

Total  Power  on  Target 

1 " Ti? 

. /d7<  |l'(r)  i2> 

rL2 

From  Eq.  (1),  we  have 

Jdr  < ! U (r ) | 2 > » V 2 | rdr  e 

Combining  these, 

U 2 a 2 

i = -5 — a_ 


(13) 


CIA) 


in  agreement  with  Eq.  (11).  Thus  the  (small-angle)  mean  irradiance  at  the 
receiver  (illuminator)  plane  is  uniform  and  independent  of  turbulence  level. 

2.  Mutual  Coherence  Function 

The  mutual  coherence  function  (MCF)  may  be  very  important  in 
analyzing  the  operation  of  a coherent  optical  adaptive  system,  and  can  be 
readily  derived  with  the  simple  assumption  of  a diffuse  target.  Ve  write 


r(pi,p?) 


(47/ 


dP,dp2  < IHpiH’MPi)  > exp 


(pj , 


Pi) 


- R2(p?.p2)]|  < exp  J\'2(Pj,P))  + <J’2*(P2»P2)] 


05) 
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where  Ri(Pi.Pi).  R2(P2»P2^  are  the  distances  from  p 
p2  respectively. 

By  the  Fresnel  approximation 


Hence 


exp 


Since  the  wave  is  incoherent  after  reflection  fron  the  diffuse 
target,  the  coherence  function  at  that  plane  can  again  be  represented  by 
the  Dirac  delta  function  as  given  in  (6).  Using  this  in  (17),  r(p;,p2) 
can  be  simplified  to 


In  the  absence  of  turbulence,  this  equation  is  entirely  identical  to  the 
Van  Cittert-Zernike  theorem  of  coherence  theory,^  which  is  identical  to  a 
result  obtained  by  Goodman  for  the  mutual  coherence  function  of  a pulsed 
optical  radar. ^ 


M.  Born  and  E.  Wolf,  Principles  of  Optics,  PerRnmun  Press.  New  York. 
1975. 

J.  W.  Goodman,  "Some  Effects  of  Target-Induced  Scintillation  on  Optical 
Radar  Performance,"  Proc.  IEEE,  53,  1688,  1965. 


To  complete  the  solution,  we  utilize  the  mean  intensity  at  the 
i .trget . For  the  focused  case,  Eq.  (10)  may  be  integrated  to  yield 


<1(P)> 


(if  !“</  1-  / r,r 


5/3 


(19) 


We  thus  have 


2 2 « 2 


r(pj,p2)  * (£)  |uol  y~J  dp  / r dr  J0  (e  pr) 


4a 

c o 


2 / 

7T  - (rj 


CXp 


^ik(pf-p^) j 


[’  T P'P  ‘ fe)  ■ 


• e 


,Uol?  / PdD  / rdr  Jo  (}.  pr)  Jo  (l  rP) 

/ 0 / 0 
? / \ ^ ^ 3 

-^•(0  -fe)  +t(^) 

o 


(20) 


trom  the  Fourier-Bessel  integral  formula. 


/-.  (£  k)  j°  (i- pr)  * • (0  % i<r'p> 


(21) 


Equation  (20)  can  then  be  simplified  and  becomes 


r(pi .P2) 

focused 


2 a 

if 


_ EL 

4a 


^ - * fef  * K (-W) 


• - - 2 

4a  2 

« <I(p)>  C 0 


fe) 


5/3 


+ 2L  (P‘'P^) 


(22) 
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Using  <I(p)>  for  the  collimated  case  in  (18)  and  simplifying 


vie  Ids 


- P2 


feW 


- 2 


fef + s M 


r(Pi.Pr)  • "Kp)*  c 
Collimated  (23) 

It  may  be  noted  that  the  MCF's  of  Eqs.  (22)  and  (23)  imply  a 
"white"  or  constant  spatial  power  spectrum  for  the  (complex)  amplitude, 
resulting  from  the  assumption  of  a diffuse  target.  The  more  interesting 
spectrum,  however,  is  that  of  the  irradiance,  as  discussed  below. 

3.  Covariance  and  Variance 

The  covariance  of  irradiance  at  the  receiver  is  given  by 


Cj  (Pi,P3> 


<I(Pl)l(P3)>  - <I(Pi)><I(P3)> 


(24) 


The  first  term,  the  correlation  function,  is  in  turn  given  by  the  extended 
Huygens-Fresnel  principle  as 


BjtPi.Pa) 


dp]dp2dp3dp4  <U(pi)i'*  (02)U(p3)U*(pi,)> 


• exp[ik(R,  - R2  + R3  - R,.)]  H(pj,p2,p3,pIi;p,,P3)  (25) 


where  H is  the  fourth  order  mutual  coherence  function  given  by 

H » <exp[C(p  i ,p  i)  + l*(02»Pl)  + ’f' <P 3 * P 3^  + ^*(p<,>P2)  (26) 

and 


R]  - 

|Pl 

" Pi 

R2  - 

Ip  l 

“ P2 

8 3 - 

1 P 3 

CO 

1 a 
1 

Ru  - 

l P 3 

- Pi, 
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After  reflection  from  the  diffuse  target,  the  fields  are  Gaussian  and 


spatially  Incoherent.  Therefore,  the  fields  at  the  target  can  be  expressed 


<U(ti)lh‘(c>2)l,(P3)bA(pe)>  “ <U(p])U-*(e2)>  <U(r>3)U'(Pi,)> 


+ <U(Pi)U'(Ps)>  <U*(f.2)U(p3)> 

“ <I(Pl)>  <I(P3)>  4(pi"P2^  ®(P3"Pi() 
+ <l(p'l)>  <1  (P 3) s «(pi"Pi,)  5(P3-P2> 


(?) 

c?y 


(27) 


We  use  this  in  Eq.  (25),  and  also  utilize  F.qs.  (A9.A10)  from  Appendix  A to 
yield 


B1(Pl.P3> 


ill 


dpjdpj  <I(pj)><I(p3)>  e 


4C  (pj-p3.P1-.--3) 


^■||dB; 


dPa  <He\)><l(.c,-})>  exp  (;  3-p  ; ) • (p;-?;)J 


. H(p  j-p  3 ,Pi_P3 ) 


(26) 


where  use  has  been  made  of 


ri-r2  - ^ [pi?  - p22  - 2(p,-p2)  pj] 

R3-R4  Zjl  |p3?  - P*.3  - 2(P3  _ P4)>] 


*ihis  assumes  that  the  rota'  po'.cr  contributing  to  the  reflected  field  at  any 
point  is  substantially  constant,  a condition  which  is  satisfied  in  the  present 
rase  of  an  extended  target  which  intercepts  the  cut  in  transmitted  beam. 
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I 


J 


We  now  combine  (8),  (24),  and  (28)  to  give  (p  = Pj-p^,  p 5 Pj-Pj): 


C_(p)  “ — — dpidp3  <I(oi)><I(p3)> 

1 h2L4  JJ 

U4C  (p.p)  \ -(ik/L)p.p 

e x -1  / + e H(p,p) 


(29) 


where  H is  given  by  Eq.  (All).  This  is  a very  important  result,  and  repre- 
sents a crucial  extension  of  the  developments  of  the  preceding  report.  In 

/ 4Cy (p.p)  S 

particular,  the  first  term  \e  A - 1/  was  not  included  in  the  earlier 
analysis  (cf.  Eq.  76,  Ref.  1)  and  represents  the  "incoherent"  contribution 
. to  the  covariance.  This  contribution  is  that  which  would  be  present  in  the 
temporally  incoherent-source  case  (Sec.  II-A4)  and  arises  from  the  log- 
amplitude  perturbation  alone.  In  contrast,  the  remaining  "coherent"  term 
(H)  involves  both  phase  and  amplitude  perturbations  multiplicatively  (Eq.  All) 
and  vanishes  in  the  incoherent-source  case. 

The  incoherent  term  will  be  seen  later  to  involve  covariance  scales 
* . 

related  to  the  "source"  size  and  /L/k,  while  the  coherent  term  basically 
involves  the  pq ("atmospheric")  and  free  space  ("target")  speckle  scales. 

Since  the  incoherent  and  coherent  terms  are  additive,  there  will  be  correla- 
tion over  the  scales  of  the  incoherent  term  even  when  the  coherent  speckles 
become  small,  as  is  observed  in  the  experiments  reported  in  a later  section. 
That  is,  the  presence  of  small,  high-contrast  speckles  will  not  obviate 
larger  correlation  scales  which  would  be  present  regardless  of  the  temporal 
coherence  of  the  source. 

We  note  that  H(p,p)  involves  both  phase  and  amplitude  perturbations, 
with  the  form  of  the  latter  strongly  dependent  on  the  strength  of  scattering 
for  a point  source  over  the  path  (Appendix  A).  We  will  accordingly  distin- 
guish "weak"  and  "saturated"  (multiple  scattering)  conditions. 

To  further  reduce  Eq.  (29),  we  utilize  Eq.  (10)  for  the  focused  case. 

The  "source"  will  hereafter  refer  to  the  spot  on  the  target. 
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and  tedious  but  straightforward  integrations.  The  result  is 

^ / n 2 V 2 f —r^/Oct  2 _ 


, .2  / “ z\z  f -r2/2ot  2 - 2(r/p  ) 

«,®  ■ f Sr  l».l*  (+)  J • ' 

o 

x | Jo  < i rp)  [ (“.«•»  .)  * „-ak'1>^  «5J)] 


The  corresponding  result  for  an  arbitrary  focus  is  obtained  by  inserting 
exp  | -r^  — j ( 1 - ) in  the  first  integral.  This  substitution 

1 

is  appropriate  in  all  similar  expressions  in  the  remainder  of  the  report, 
so  that  the  discussion  can  be  confined  to  the  focused  beam  case  without 
loss  of  generality. 

These  results  (30)  are  five-fold  integrals  but  not  formidable;  their 
physical  interpretation  and  numerical  evaluation  will  be  discussed  further 
in  later  sections. 

The  variance,  normalized  by  the  square  of  the  mean  irradiance,  is 
given  by  Eqs . (30)  and  (A12): 

r rn\  . -00  _ ~ 2 / 2 „ -W _ / ~ \5/3 


CIW  k2  f 

*n  U 4a  “/L4  L2  J 

0 0 0 


-r2/2a  2 - 2(r/p  )" 
o o 


C k 

r 4cv(c)  i 

j odo  Jo  ( L rp) 

L 2e  x - 1 J 

In  order  to  identify  the  incoherent  and  coherent  contributions  respectively, 
the  final  bracket  may  be  written 
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It  is  of  interest  to  compare  this  variance  result  with  that  obtained 
in  the  preceding  report  under  the  assumption  that  phase  perturbations  domi- 
nate. In  that  case,  the  normalized  variance  was  found  to  be  identically 
unity.  From  the  present  complete  analysis,  we  see  that  the  phase  perturba- 
tions do  not  explicitly  enter  the  variance  expression,  the  phase  having  been 
totally  randomized  by  the  diffuse  target.  In  the  limit  of  weak  turbulence 
(C^  = 0),  the  integrations  in  Eq.  (31)  will  yield  a variance  of  unity,  physi- 
cally owing  to  the  diffuse  target  (free-space  speckles);  in  the  limit  of 
strong  turbulence,  it  will  also  be  unity  (C^  = 0 for  P>PQ  c e,  where  e is 
vanishingly  small).  For  intermediate  values  of  turbulence,  the  variance  will 
be  somewhat  greater  than  unity,  as  a result  of  non-negl igible  amplitude  per- 
turbations (C  ).  These  physical  interpretations  will  be  discussed  in  detail 
in  a later  section. 

We  may  often  assume  in  Eqs.  (30,31)  that  4C^<<1,  so  that  (2e  -1) 

% 1 + 8C^.  This  is  a good  approximation  in  weak  turbulence  and  a rough  but 
useful  approximation  in  strong  turbulence,  and  may  substantially  simplify 
further  analysis.  For  example,  Eq.  (31)  with  (A4)  becomes 


These  expressions  clearly  show  the  increment  from  unity  variance  related  to 
the  log  amplitude  perturbation.  Large  values  of  and  p^,  corresponding 
to  small  target  spot-size,  result  in  greater  variance  (stronger  scintilla- 
t ions) . 


4.  Incoherent  Case 


In  order  to  derive  the  covariance  for  a temporally  incoherent 
source,  we  return  to  Eq.  (25)  and  write 

Bj  (P1.P3)  - Wf  dp  jdp2dP3dc>4  5SSS  dkidk2dk3dki, 

* i (Pi  (p2»^2)^3(P3»*c3)lJi»(P4t*ci))> 

• exp  { i (kiRi-k2R2+k3H3-kt,R4)^  • H (pj-  • •Pi'iPi.Ps)  (34) 

We  then  note  that  the  incoherence  of  the  source  implies  that  the  integrals 
over  kj_4  will  yield  zero  except  for  (kj.pj)  = (k2,P2)  and  (k3>  o3)  = (k^.p^) 
This  drops  the  phase  perturbation  out  of  H,  by  virtue  of  Eq.  (A9),  and  we 
are  left  with 


_ _ _ _ _ ^*C  (p  -p  ,p) 

Bj(p)  dptdo3  <I(Pl)><I(p3)>  • e 


(35) 


The  corresponding  covariance  is  then  recognized  as  the  "incoherent  part"  of 
the  general  covariance  of  Eq.  (2S). 

We  also  point  out  that  this  result  is  easily  shown  to  be  consistent 
with  the  results  of  Ref.  7,  in  which  a spatial-frequency  domain  treatment  is 
given  for  a target  under  ambient  illumination. 

5.  Receiver  Smoothing 

In  this  section  we  treat  the  smoothing  of  diffuse-target  scintilla 
tions  by  a finite  receiver  aperture. 

We  define  the  receiver  weighting  function  as  W(p,D): 


W(p,D) 


1 ^ Ip  I 1 § 

0 if  \p\  > | 


(36) 


7.  S.  F.  Clifford,  G.  R.  Ochs,  and  Ting-i  Wang,  "Optical  Wind  Sensing  bv 
Observing  the  Scintillations  of  a Random  Scene",  Applied  Optics,  14, 
2844,  December  1975. 


The  mean,  smoothed  (aperture-integrated)  receiver  intensity  is  simply  (Eq.ll) 


<Is>  = J dp  W(p,D)<l(p)> 

- ^ u 2 *4 

4 o l2 


The  mean-square  intensity  is 

<xs2>  " SS  dpjdp2  <I(pi)  l(p2)>W(p;,D)W(p2,D) 

and  the  variance  is 

°j2  “JSdpjdpz  Vs’(pi,D)K(p2,D)  C (pltp2) 

s 

We  let  Pi_P2=p,  pi+p2=2p'  and  integrate  over  p': 

M(p,D)  £ J dp'W(p’+p/2,D)K(p,-p/2,D) 

- D2/2  jcos'1  (p/I))-  (p/D)  [l-p2/D2]  if  p < D 
= 0 if  p > D 
Letting  x = p/D,  we  have 


(37) 


(38) 


(39) 


irD^  f 

"4  3 


dx  M(x)  Cj  (Dx) 


lxl<l 


With  M(x)  = 2/tt  [cos  * x - x /l-x2  3 


(40) 


and  Cj  is  given  by  Eq. 


(30). 


In  order  to  obtain  Lhc  finite-receiver  covariance,  we  define  variables 
Tj.r^  as  in  Figure  2.  Then 


<1  > - J <I(Pi+  rj  )>W(r  j ) dr  j 

“ X <l(P2+P+ri)>W(rI)dri 
<IS2>  " X <I(P2+r2>>W(72)d72 
and  finally 


Cj  (p)  = irDV4  J dx  M(x)  Cx  (p  + Dx) 


I 


6.  Higher  Moments  and  Probability  Distribution  of  Irradiance 

The  phase  and  amplitude  perturbations  imply  exponential  and  log 
normal  irradiance  statistics  respectively.  It  is  of  great  interest  to  define 
the  composite  probability  distribution  as  a function  of  the  pertinent  param- 
eters. Analytical  and  qualitative  considerations  will  be  reviewed  in  this 
section;  a complete  treatment  is  quite  complex  and  has  not  yet  been  accom- 
plished . 

We  approach  the  probability  distribution  through  its  moments.  From 
Eq.  (28),  we  write  the  second  moment  of  irradiance  as 


<I2(p)>  = 2/tt2L4JJ  dPldp3  <l(pi)><l(p3)>  e 


4CV< ioi~P3  f > 


This  may  be  generalized  to  the  nth  moment: 


l"(p)>  -•  n!  (^2-)  fdpi...dpT„r  n <I(p-.)>]  exp  [2  T C 1 

J J L i=l  J L ijfj  xijJ 


We  write  this  as 


where 


Note  tnat  c 


The  interpretation  of  these  expressions  is  facilitated  if 
then  simplifies  to 

4jJdp1d03<I(F1)><i(p3)>  Cv(!p,-r,j) 


where  the  denominator  is 


the  integral  in  the  numerator 


MluU 


For  saturated  scintillations,  Eqs.  (10)  and  (AS)  give 


4 ««■>  f *• | - 1 - * ftJW’-fe  - 1 1 w 5/3  <>”> 

J * o o a o 


We  then  have  the  second  moment  from  Eq.  (45): 


<I2(p)>  "2<I(p)>2  F2 


and  the  nth  moment  is 


<tn>  % n!  <!>*'  \ 1 + (F2-l)| 


Hence  nonunity  F.  represents  the  departure  from  an  exponential  intensity 
distribution.  Note  that  F2  will  be  unity  for  weak  turbulence  (C  *v  0),  or  for 
a large  target  spot  resulting  from  a small  transmitter  (small  a q) , deliberate 
defocusing,  or  strong  turbulence  (small  pQ). 

Eq.  (52)  is  only  a good  approximation  for  n such  that  nC  «1.  It  is 
tempting  to  derive  the  probability  distribution  using  (52)  with 

<e~sI>  - l-s<I>  + + • • • (53) 


Defining  <I>  « 1^,  the  result  is 


p,d)  - r . 
o 


|l  + (F2-l)  £l  - ■ip-  + | + (higher-order  terms) 


where  the  "higher-order  terms"  recognize  the  limitations  on  (52) . Although 
this  result  appears  to  give  a first-order  correction  to  the  exponential  dis- 
tribution, it  is  apparent  that  the  higher  moments  play  a crucial  role  in  the 
probability  distribution  since  this  apparent  correction  does  not  agree  with 
physical  reasonings. 

Specifically,  let  us  take  the  heuristic  viewpoint  that  the  irradiance, 
at  least  for  small  C^,  is  comprised  of  the  sum  of  two  random  processes,  with 
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exponential  and  log  normal  distributions  respectively.  The  composite  distri- 
bution is  given  by  the  convolution  of  these  distributions,  with  the  qualita- 
tive effect  shown  in  Figure  3.  This  prediction  agrees  with  experimental 
results  discussed  in  Section  IV,  but  not  with  the  expression  in  (34). 

Work  is  continuing  on  the  derivation  of  analytic  expressions  for  the 
probability  of  irradiance. 


7.  Physical  Interpretation  of  Variance  and  Covariance 


The  physical  interpretation  of  the  variance  and  covariance  as 
expressed  in  Eqs.  (30)  and  (31)  respectively  reveals  the  mechanisms  Involved 
in  speckle  propagation  through  turbulence  in  a particularly  illuminating 


The  variance,  Eq.  (31),  may  be  expressed  as 


which  in  many  cases  can  be  written 


'v  l + 8y  o 2 
X 


In  these  expressions,  a2  is  the  log  amplitude  variance  for  a point  source 
propagating  over  the  target-receiver  path,  and  y (<  I)  is  the  "source  smooth- 


ing factor”  associated  with  the  nonzero  spot  size  (9)  on  the  target.  Mathe- 


matically, y relates  to  the  integration  over  (p)  in  Eq.  (31),  with  the  spot 


size  (S  “ or  ) coupled  through  r.  For  a small  transmitter  (or 
defocused  spot),  or  strong  turbulence,  the  target  spot  will  be  S * /L/k  and 


Y will  be  less  than  unity.  In  the  limit  of  strong  turbulence,  y ^ 0 while 
for  weak  turbulence,  a2  ” 0;  either  extreme  results  in  a unity  variance. 


for  weak  turbulence,  o2  » 0;  either  extreme  results  in  a unity  variance. 

A quantitative  estimate  of  y may  be  obtained  as  follows.  From  pub- 

8 . . . . L 


liahed  results  for  the  variance  with  an  extended  Incoherent  source,  we  know 
that  in  weak  scattering 


8.  R.  F.  Lutorairskl,  et  al.  Degradation  of  laser  Systems  by  Atmospheric 
Turbulence,  HAND  Report  R-1171-ARPA/RC,  June  1973. 
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Log  Normal  Distribution 
(Amplitude  Perturbations) 


Exponential  Distribution 
(Phase  Perturbations) 


Composite  Distribution 


Figure  3.  Qualitative  Probability  Distribution  of  Irradianre 


(56a) 


The  smoothing  in  the  strong  scattering  case  is  conjectured  to  be 
-2 


( -) 
v Pc  ' 


(56b) 


This  source  or  target-spot  smoothing  is  analogous  to  finite-receiver  smooth- 

. 9 

ing. 

The  expressions  corresponding  to  (55)  for  an  incoherent  source  are 
(Sec.  11A4) : 


4y  o ‘ 

2 = o X 


I 


- 1 


(57) 


y ^ . y 

cl*  ^ Ay  o * 
X X 


(57a) 


It  may  been  seen  that  the  unity  baseline  and  half  the  remaining  term  in  (55a) 
arise  from  the  source  coherence  or  speckle  phenomenon,  whether  controlled  by 
the  target  (weak  turbulence)  or  atmosphere  (strong  turbulence).  Alternative- 
ly,  we  may  observe  that  the  unity  term  arises  strictly  from  phase  randomiza- 
tion, with  associated  exponential  irradiance  statistics  (Ref.  1),  while  the 
remaining  term  represents  the  effects  of  nonzero  amplitude  perturbations, 
with  associated  log  normal  statistics.  This  supports  the  heuristic  probabil- 
ity arguments  of  the  previous  section. 

To  further  illustrate  these  considerations,  let  us  consider  the  large, 
focused  transmitter  condition  (a ^>>\/h/k)  (Figure  4).  We  first  consider  the 


weak  turbulence  condition  (p0>aQ>/l./k) , in  which  the  transmitter  is  able  to 

operate  in  an  essentially  diffraction-limited  manner.  The  target-spot  is 

small,  i.e.  such  that  the  receiver  is  in  the  far-field  of  the  spot,  and  the 

situation  is  much  like  that  of  a point  source  (o2<<l,y=l).  There  is  however 

X 

one  major  difference  from  the  point-source  case:  we  assume  a good  statistical 


9.  R.  S.  Lawrence  and  J.  W.  Strohbehn,  "A  Survey  of  Clean-Air  Propagation 
Effects  Relevant  to  Optical  Communications",  Proc.  IEEE,  5(5,  1523, 
October  1970. 
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Turbulence  Strength  (Vpq) 
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o o 

o 2 % 0.6 

X 

y < l 


{a  >/L/k>p  } 
o o 
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X 
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Figure  4.  Normalized  Variance  of  Trradiance  for  Increasing  Turbulence 


(ensemble)  average  of  the  field,  which  as  a practical  matter  means  that  the 

target  is  moved  slowly  so  that  the  (large)  speckles  sweep  past  the  receiver. 

(If  the  transmitter,  target,  and  receiver  were  rigidly  fixed,  the  receiver 

would  see  part  of  one  speckle,  and  the  turbulence  fluctuations  would  be  that 

of  a point  source,  with  purely  log  normal  statistics  and  a variance  of 
Aov2  <v  , 

e *■  - 1 'v  4 a *).  The  resultant  variance  is  essentially  unity,  with  small 
turbulence-caused  fluctuations  superimposed  on  the  basic  speckles  that  would 
be  present  in  the  absence  of  the  atmosphere.  The  irradiance  statistics  will 
be  nearly  exponential. 

We  now  let  the  turbulence  strength  grow  until  "saturation”  or  multiple 
scattering  sets  in.  The  peak  value^  of  a 2 is  known  to  be  approximately  0.6, 


occurring  when  is  on  the  order  of  /L/k.  The  target  spot  will  no  longer 
be  diffraction  limited  (a  >p  "'I  /L/k)  ,and  hence  y will  be  less  than  unity.  The 
exact  v»li!e  of  the  peak  variance  can  only  be  determined  from  a full  numerical 
evaluation  of  F.q.  (31),  including  the  detailed  effects  of  saturation  on 

(p);  the  statistics  will  be  a composite  of  log  normal  and  exponential  dis- 
tributions. As  the  turbulence  strength  grows  further,  o 2 decreases  ("super- 
saturation")  to  the  licit  0.17;  also  the  target-spot  grows  indefinitely,  so 
that  y decreases  to  2ero  in  the  limit.  The  variance  therefore  decreases  back 
to  unity  (Fig.  4),  with  exponential  statistics. 

We  now  consider  the  covariance  case,  again  for  a focused  transmitter. 

We  refer  again  to  Eq.  (30).  For  the  coherent  (H)  term,  if  is 

neglected  and  either  a or  o is  much  smaller  than  the  other  and  also  smaller 
o o 


than  /L/k,  the  corresponding  limitation  on  r in  the  integrand  couples  throjgh 
the  Bessel  term  and  exponential  to  imply  a similar  limitation  on  p.  Referring 
to  Eq.  (All),  the  p-terms  cancel  and  H = e ' . This  gives  a closed- 

form  covariance  result  which  is  that  obtained  with  the  assumption  of  joint 
gaussian  field  statistics  (Ref.  1): 


(p) 

n 


-p2/2a  2 - 4(p/p  )5^3 


Hence  the  covariance  arising  from  the  coherent  term  in  Eq.  (30)  is  generally 
exponential  with  a scale  equal  to  the  smaller  of  a0  or  p0  (target  and  atraos- 
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pheric  speckles  respectively).  In  the  event  that  a0  and  pQ  arc  comparable, 
the  complete  expression  must  be  used  and  the  field  statistics  are  marginally 
(but  not  joint)  gaussian.  In  the  event  that  is  non-negl igible , the 
covariance  scale  relating  to  the  coherent  term  is  largely  unaffected. 

We  now  consider  the  incoherent  term  in  (30) . If  the  target  spot  is 
small  compared  to  /L/k,  i.e.  a0  and  Po>>/L/k,  the  large  range  of  r in  the 
integrand  will  couple  through  to  hold  p small  in  C^.  Hence  the  covariance 
scale  will  be  determined  by  C (p,0)  and  will  be  /L/k  as  expected  for  a small 
target  spot.  However,  for  a large  target  spot  (S),  the  p range  in  the  inte- 
grand will  be  large  and  the  p-scale  in  will  be  "smeared"  to  the  order  of  S. 
This  effect  is  consistent  with  the  incoherent-source  results  of  Kef.  7,  and 
can  be  argued  from  simple  geometry:  for  two  receivers  separated  by  a small 

distance  relative  to  the  source  size,  the  source-turbulence  configuration  over 
the  paths  to  each  respective  receiver  will  be  highly  correlated.  The  statis- 
tics associated  with  the  incoherent  term  are  log  normal,  although  the  log 
amplitude  is  not  joint  normal  except  for  weak  fluctuations^ 

We  now  summarize  the  above  reasonings  with  respect  to  particular  param- 
eter realms: 


a)/l7k<ao,p0 


In  this  case,  the  target  spot  is  small,  the  receiver  is  in  the  far 
field  of  the  spot,  and  y ” 1.  The  covariance  contains  both  a phase  and  ampli- 
tude term  as  shown  in  Figure  5a.  The  amplitude  term  has  a variance  contribu- 
4ov^ 

tion  of  2(e  -1)^8  o 2 and  a scale  of  /L/k,  while  the  phase  term  has  a vari- 

ance of  unity  and  scale  equal  to  the  smaller  of  (ao,PQ).  The  latter  ',cale 
represents  the  speckle  size,  as  determined  by  the  target  spot  size  (L/ka0  or 
L/kpo).  Region  II  is  dominated  by  free-space  speckles,  such  that  the  associ- 
ated temporal  spectrum  is  determined  by  target-receiver  motion  and  the  sta- 
tistics* relate  to  the  exponential  irradiance  distributior. . In  Region  1,  the 


*These  statements  regarding  statistics  related  to  covariance  realms  do  not 
necessarily  imply  joint  log  normal  or  joint  exponential  statistics  for  two 
points  separated  by  p,  but  rather  relate  to  the  single-point  statistics  of 
the  corresponding  mechanisms. 

10.  J.  R.  Kerr,  et  al.  Propagation  of  Mult ivavelength  Laser  Radiation 

Through  Atmospheric  Turbulence,  RADC-Tf -76-49,  March  1976,  (A023202) . 
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increment  from  unity  represents  the  amplitude  or  log-normal-like  contribu- 

* 

tion,  with  a temporal  spectrum  related  to  the  transverse  wind.  Since  the 
fall-off  (Region  II)  of  the  phase  term  is  basically  exponential,  while  the 
amplitude  term  falls  off  more  slowly  and  with  an  oscillatory  nature,  there 
will  exist  a tail  (Region  III)  relating  entirely  to  the  latter. 

b)a  </L/k,o 
o u 


In  this  case,  the  target  spot  is  large,  the  receiver  is  in  the  near 
field  of  the  spot,  and  y<l.  The  amplitude  term  has  a contribution  of  8yo^2 , 
and  a scale  which  is  caused  by  geometrical  considerations  to  be  on  the  order 

of  the  target-spot  size  (L/ka  >/L/k) . The  phase  term  starts  at  unity  and  has 

0 

a scale  of  a^,  representing  the  speckle  size.  Region  11  is  the  log-normal- 
like  region,  with  a temporal  spectrum  related  to  the  transverse  wind;  Region  1 
includes  the  contribution  related  to  the  exponential  radiance  distribution, 
with  temporal  behavior  related  to  target-receiver  motion. 


c)p  <>L/k,a 
o o 


In  this  case,  the  target  spot  is  large  owing  to  atmospheric  degradation 

of  the  transmitter  wavefront,  the  receiver  is  in  the  near  field  of  the  spot, 

and  y<l.  The  considerations  are  similar  to  those  for  the  preceding  case, 

except  that  cq  replaces  a^  in  determining  speckle  size.  Since  the  phase  tern 

is  now  dominated  by  this  turbulence-related  quantity,  it  also  has  a temporal 

behavior  related  to  the  transverse  wind. 

A useful  way  to  summarize  the  above  considerations  is  in  tabular  form 

as  in  Ref.  1.  The  further  developments  of  the  theory  contained  herein  require 

a modification  of  Tables  III  and  IV  of  that  report,  as  follows. 

For  the  focused  case,  the  independent  parameters  of  interest  remain 

a , •'TTk,  and  p . Since  there  will  now  be  two  additive  terms  in  the  covari- 
o o 

ance,  there  will  in  general  be  two  covariance  scales  (Figure  5);  the  scale 
of  the  coherent  (phase)  term  will  dominate,  owing  to  the  weaker  fluctuations 
associated  with  the  incoherent  term,  but  in  cases  (b,c)  of  Figure  5 there 
will  be  a significant  correlation  "tail"  beyond  these  scales. 
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The  results  are  given  in  Table  1.  In  parameter  realms  A and  B,  in 
which  the  receiver  is  in  the  far  field  of  the  small  target  spot,  we  distin- 
guish two  cases.  In  case  1,  we  assume  strict  mechanical  stability  between 
the  transmitter,  target,  and  receiver,  such  that  the  (very  large)  speckles 
are  not  moving.  Each  receiver  will  then  sample  an  unvarying,  snail  portion 

of  a speckle  and  the  only  fluctuations  will  be  those  arising  as  from  a point 

9 

source  in  turbulence;  this  well-known  situation  involves  a different  defi- 
nition of  statistical  averaging  than  that  which  is  implied  in  the  foregoing 
analysis.  In  case  2,  the  ergodic  ensemble  average  is  used  which  corresponds 
to  movement  of  the  speckles  past  the  receiver.  As  a practical  matter,  an 
imperfectly  "fixed"  target-transceiver  configuration  will  behave  somewhere 
between  these  two  extremes;  rotation  or  overt  motion  of  the  target  however 
results  in  the  unambiguous  realization  of  case-2  statistics. 

For  the  unfocused  case,  with  an  arbitrary  spot  size  (S)  assumed  to  be 
>/L/k,  we  have  three  corresponding  realms  as  shown  in  Table  II.  The  special 
case  of  collimation  is  given  by  S = aQ. 

Finally,  we  consider  physically  the  covariance  characteristics  for 
finite  receivers.  The  result  expressed  in  Eq.  (42)  infplies  that  the  covari- 
ance scale  for  both  the  coherent  or  incoherent  terms  in  C^(p)  will  be 
smeared  to  the  order  of  the  receiver  size,  with  a falloff  whose  abruptness 
depends  on  the  point-receiver  covariance  scales  relative  to  this  size.  This 
behavior  is  obvious:  separations  less  than  the  receiver  size  imply  overlap- 

ping receivers.  The  variance  for  a finite  receiver  will  be  smoothed  in  ac- 
cordance with  the  effective  number  of  independently  scintillating  patches  in 

g “7/3 

its  aperture;  the  dependence  for  the  incoherent  term  is  (D/>T./k)  in  the 

absence  of  multiple  scattering  and  a function  of  (D/p^)  for  strong  scattering; 

the  dependence  for  the  coherent  term  has  not  been  explicitly  obtained  but  may 

“ 2 

be  expected  to  be  approximately  ^ D/(a ^ or  pQ)^j 

(smaller) 

We  point  out  that,  for  a receiver  large  enough  to  smooth  the  speckles 
represented  by  the  coherent  term,  these  speckles  become  simply  a residual 
noise  source  and  the  behavior  is  otherwise  as  if  the  source  were  incoherent. 
Referring  to  Figure  5,  the  tail  of  the  covariance  term  then  has  a more  sig- 
nificant effect  than  the  smaller  scales. 
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8.  Time-Lagged  Covariar.ee 


The  spatial  spectrum  of  irradiance  may  be  found  by  taking  the 
Fourier  transform  of  the  covariance  (Eq.  30).  The  temporal  spectrum  requires 
the  time-lagged  covariance  function,  which  may  be  obtained  by  repeating  the 
procedure  of  Sec.  1I-A3  incorporating  the  time-delay  (t)  between  irradiance 
measurements  at  p^  and  p^  respectively.  We  have  the  following,  corresponding 
to  Eqs.  (24,  25,  and  27): 

Cj(Pl  .P3.T)  - <I(pj)I(p3,t)>  - <I(pj)><1(p3;t)>  (59) 

BI<Pl»P3»T)  “ (k/2i:L)4  dp1dp2dp3dp1,  <l'(Pi)U  (p2)U(p3,t)U  (Pi,,t) 

exp  { ik(R1-R2+R3-R,)i  H(r  j , P 2 ,P3 .Pi, Ipi  ,p2 ; t)  (60) 

* _____  2 
<U(0!)-  • -U  ( p t,  ,t)>  = <I(p2)><I(p3,T)>c(p1-p2)6(p3-.-^X4-/k2) 

+ <U(pj)U  (cj,t)><U  (p3)U(c3,T)>i(.,1-Pi,)6(p2-p3)  (4r/k2)2  (61) 

where  <I(p1>t)>  = <I(ci)> 

The  problem  is  now  complicated,  relative  to  the  previous  (t  = 0)  case, 
by  the  time-delay  term  in  the  second  part  of  Eq.  (61).  That  is,  in  order  to 
insert  Eq.  (61)  into  (60),  we  must  return  to  Eq.  (3)  for  E’(p)-  This  in  turn 
introduces  the  time-delayed  wave  structure  function  D^(r,r).  The  ensuing 
integrations  will  again  be  omitted;  the  result,  which  utilizes  the  results  of 
Appendix  B and  corresponds  to  Eq.  (30),  is  ' 


C (p.T)  = C + C 

i ij  i2 


C_  = 2/it(k2/L6)  |U  |4  (a  2/2): 

I J O 0 


oo 

j"  rdi 


-r2/2a  2 - 2(r/p  )^^ 
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Work  is  currently  underway  on  the  interpretation  of  these  expressions. 
In  particular  the  value  of  -C  /it  at  t = 0 will  give  insight  into  the  incre- 
mental motion  vs.  evolution  of  the  scintillation  pattern. 

II-B.  Complex  Targets  and  Sources 

The  preceding  considerations  were  based  on  the  assumption  of  a perfect- 
ly diffuse  target.  In  order  to  predict  the  effects  of  real  targets,  such 
features  as  glints,  boundaries,  and  phase  correlation  areas  (specular  compo- 
nents) must  be  included.  In  addition,  the  imperfect  temporal  and  spatial 
characteristics  of  real  lasers  may  have  an  influence  and  should  be  understood. 

These  topics  are  currently  under  study,  and  introductory  considerations 
will  be  briefly  presented  in  this  section.  Ultimately , it  will  be  fruitful 
to  extend  this  work  to  more  general  target  models  and  scenarios,  including 
particular  wavelengths  of  interest,  moving  target/transceiver  configurations, 
and  slant  paths;  and  to  deal  directly  with  turbulence-  and  target-driven 
effects  on  coherent  adaptive  transmitter  svjeems. 

1.  Target  Truncation 

The  simplest  example  of  a nonideal  target  is  a diffuse  surface  with 
boundaries.  To  the  extent  that  the  boundaries  (which  might  represent  sub- 
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features  of  a more  complex  target)  appreciably  truncate  the  target-spot,  there 
will  be  a component  of  high  spatial  frequencies  at  the  target  and  low  frequen- 
cies at  the  receiver,  l.e.,  there  will  be  larger  speckles  at  the  receiver, 
which  may  represent  a significant  noise  source  in  a finite  aperture. 

We  again  treat  a focused  beam,  where  the  extension  to  a more  general 
focal  condition  is  straightforward  as  in  the  preceding  developments.  We 
assume  a round  target  of  radius  (a).  Referring  to  Eqs.  (8),  (19)  we  have  for 
the  mean  receiver  irradiance: 

k2U  2a  2 f 3 f -r2/4ci  2 -(r/p  )5^3 

<HP)>  - * P-  l Pdpl  rdr  J (k  pr/JJe  0 0 (64) 


For  a more  general  target  shape,  (S),  we  can  write 

- k2U  a t “ -r2/ 4a  2 - (r/i 

<I(p)>  “ °-  -9-  l rw(r)  e 


w(r)  = (a  2 / 2 ) ( 

kar/L 

Similarly,  we  write  the  covariance  from  Eq.  (29)  with  appropriate 
limits  on  Pi,3,  and  let  pj-p3=p,  pj+p3=  2R: 


(k/L)2  U 


The  result  is  somewhat  more  complex  than  that  of  Fq.  (30)  for  an 
infinite  target.  A full  interpretation  and  numerical  evaluation  remains  to 
be  completed.  However,  the  physical  implications  can  be  generally  deduced. 

Let  us  assume  that  (a)  is  much  smaller  than  the  spot  size  determined  by  a : 

o 

Coherent  term:  the  speckles  scale  arising  from  the  target  will  be 

'vL/ka>>aQ.  Phase  perturbations  will  cause  a normalized  variance  of  unity; 
amplitude  perturbations  may  have  a stronger  additional  effect  than  in  the 
infinite-target  case,  since  the  total  power  hitting  the  target  may  now  vary. 

Incoherent  term:  The  limited  size  of  the  target  may  reduce  or  elimi- 

nate spot-size  smoothing  of  the  receiver  scintillations,  which  will  be  further 
enhanced  by  variations  in  the  total  power  hitting  the  target. 

We  note  here  a very  interesting  prediction  that  may  be  subjected  to 
numerical  and  experimental  verification,  and  which  could  have  profound  effects 
on  optical  system  operation.  Suppose  that  the  target  (or  some  feature  or.  the 

target)  is  smaller  than  the  scales  of  interest,  i.e.,  a<P  ,>'T7k.  We  may  then 

0 

utilize  <5-functions  on  p,R  in  Eqs.  (68,69)  and  find  that  the  variance  is  given 
by 


n 


(70) 


where  we  have  again  deliberately  separated  the  contributions  from  the  incoher- 
ent and  coherent  terms  respectively.  This  expression  is  consistent  with  Eq . 
(55)  with  Y”I , as  expected. 

We  now  note  that  for  a small  source  in  strong  turbulence  (c  <*'u/k), 

11  ° 

the  saturated  variance  for  the  incoherent  term  is  unity.  This  implies  that 
4o  2 

e X « 2.  Hence  under  the  present  condition,  wherein  the  effective  target 
spot  is  constrained  to  be  small  through  target  truncation,  strong  turbulence 
may  result  in  a total  normalized  variance  of  3,  which  is  indeed  substantial. 
This  effectively  comes  about  through  a combination  of  varying  phase  perturba- 
tions (speckles)  and  saturated  amplitude  scintillations  from  the  target-spot 
to  the  receiver.  (It  is  of  course  assumed  that  the  diffuse  target  has  suffi- 
cient extent  to  create  speckles,  which  are  modulated  by  the  turbulence  and/or 
target  motion.)  The  variance  of  3 can  in  fact  be  simply  related  to  the  fol- 

11.  R,  JTinte,  "Electromagnetic  Beam  Propagation  in  Turbulent  Media",  Proc. 
IEEE,  63,  1669,  Dec.  1975. 
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lowing: 

Let  random  independent  processes  (x,y)  have  normalized  variances  of 

unity. 

Let  z = xy 

Then 

°z  _ <x2y2>  _ j _ <x2xy2>  _ j 

<z>2  <xy>2  <x>2<y>2 

“ 3 


If  we  identify  (x,y)  with  (saturated)  amplitude  and  phase-induced  irradiance 

fluctuations  respectively,  we  derive  the  expected  total  variance  of  3. 

Further  insight  can  be  gained  by  considering  the  case  of  a one-way  path 

from  a finite,  diffuse  source.  We  assume  that  this  source  has  a diameter  D 

and  uniform  irradiance  1 . From  Eq.  (29)  and  the  overlap  integral  such  as 

o 

used  in  Sec.  II-A5,  we  immediately  write  the  covariance  as 


C^P) 


ik(D/l)  p-x 


(' 

2,  j 

( 4Cx(Dx,p)  ) 

1 xdx 

X de  M(x)  1 

{ \ e - 1 f 

Dx,p)  S 


+ e " H(Dx, 

where  x = p/D,  6 is  the  angle  between  p and  p,  and 

] 


M(x)  » 2/tj  [ 
Since  <I(p)>  » I 


cos  * x - x/l-x2 


(71) 


(72) 


n D2 


from  Lambert’s  Law,  the  normalized  variance  is 


uL2 


C r 4Cx(Dx)  I 

I xM(x)  L2e  - 1 J 


dx 


(73) 
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Kor  D less  than  the  characteristic  scale  of  C^,  Eq.  (73)  reduces  to  Eq.  (70). 

It  is  at  first  perhaps  curious  that  the  same  result  for  the  variance 
was  obtained  for  both  the  small  target  and  the  small  diffuse-source  cases. 

This  suggests  that  any  fluctuation  (scintillation)  of  the  total  energy  on  the 
target  is  not  taken  into  account  in  Eq.  (68).  This  is  indeed  the  case;  as 
discussed  in  conjunction  with  Eq.  (27),  the  assumption  of  gaussian  field  sta- 
tistics after  reflection  from  the  target  implicitly  drops  such  fluctuations. 

The  small-target,  strong-scintillation  situation  represents  a practical  case 
where  this  assumption  may  not  be  valid  and  a more  complete  treatment,  includ- 
ing irradiance  statistics  on  the  target,  should  be  used.  This  is  turn  implies 
that  the  receiver  variance  may  be  greater  even  than  3. 

2.  Glints 

An  analysis  of  the  mean  irradiance  pattern  at  the  receiver  for  an 
arbitrary  number  of  glints  on  the  diffuse  target  was  given  in  the  preceding 
report.  We  have  subsequently  derived  the  mutual  coherence  function  and  co- 
variance  under  the  assumption  of  joint  gaussian  field  statistics.  The  limi- 
tations of  this  assumption  have  been  pointed  out  in  the  discussions  of  this 
report,  and  an  extension  of  the  theory  is  underway. 

3.  Nonideal  Sources 

The  temporal  and  spatial  coherence  characteristics  of  real  lasers, 
especially  in  the  case  of  ma  y high  power,  infrared  units  of  interest,  cannot 
be  neglected.  The  nonzero  spectral  width  covered  by  one  or  a series  of  laser 
lines  (transitions)  may  result  in  sharply  reduced  speckle  contrast  in  the  co- 
herent term  of  the  variance,  and  possibly  even  spectral  smoothing  of  the  inco- 
herent term  in  the  presence  of  strong  turbulence.  These  effects  will  be 
observed  if,  as  is  the  case  in  the  initial  experiments  described  in  Section  TV, 
the  spectral  width  leads  to  a "coherence  length"  which  is  less  than  the  path- 
length-dif ferences  involved  in  the  optical  system-target-atmosphere  geometry. 

In  addition,  transverse  mode  structure  can  increase  the  focused  target-spot- 
size  and  therefore  affect  the  received  field. 

The  effects  of  finite  spectral  width  on  speckles  in  the  absence  of 
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turbulence  have  been  heavily  researched  in  the  past  few  years.  The 

extension  of  the  current  theory  to  include  both  temporal  and  spatial  coher- 
ence effects  in  turbulence  is  currently  underway. 

A heuristic  approach  is  suggested  as  follows.  Ke  separate  the  coherent 
and  incoherent  parts  of  the  variance  (Eq.  55): 


n 


A*2°x2  ( ^2°  2 \ 

Y i e + \e  x - 1 / 


(74) 


where 

Yj  “ reduction  in  speckle  contrast  owing  to  source  spectral  width 
Yj  “ source  (target-spot)  smoothing 


While  it  is  recognized  that  Yj  will  vary  with  the  details  of  transmitter- 
target-receiver  geometry,  including  spot-size,  this  simple  approach  is 
shown  to  give  reasonable  qualitative  agreement  with  experiments  in  Section 
IV. 

III.  Numerical  Evaluations 

The  results  of  Section  11  for  such  statistical  descriptors  as  the 
covariance  of  irradiance  are  generally  in  the  form  of  multiple-dimensioned 
integrals.  While  amenable  to  meaningful  physical  interpretation,  accurate 
quantitative  predictions  from  these  expressions  requires  numerical  evaluation. 
The  purpose  of  this  section  is  to  present  a simplifying  technique  which  has 
made  the  problem  tractable  by  greatly  reducing  the  time  required  on  a com- 
puter of  reasonable  capacity  and  speed.  Preliminary  results  are  given  and 
compared  with  experiments  in  Section  IV. 


12.  J.  C.  Dainty,  "Some  Statistical  Properties  of  Random  Speckle  Patterns  in 
Coherent  and  Partially  Coherent  Illumination'',  Optica  Acta  l_7,  761,  1970. 

13.  C.  Parry,  "Sou.e  Effects  of  Temporal  Coherence  on  tac  First  Order  Sta- 
tistics of  Speckle",  Optica  Acta  2\_,  763,  1974. 

14.  H.  M.  Pederson,  "On  the  Contrast  of  Polychromatic  Speckle  Patterns  and 
its  Dependence  on  Surface  Roughness",  Optica  Acta  22,  15,  1975. 

15.  H.  M.  Pederson,  "Second-Order  Statistics  of  Light  Diffracted  from  Oaus- 
sian.  Rough  Surfaces  with  Applications  to  the  Roughness  Dependence  of 
Speckles",  Optica  Acta  22,  523,  1975. 
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III-A.  Basic  Technique 

The  five-fold  integrations  required  in  expressions  such  as  (30) 
initially  appear  rather  formidable.  The  use  of  Monte  Carlo  techniques  is 
hampered  by  the  form  of  the  integrand  and  the  interdependency  of  the  inte- 
gration limits.  An  analytic  technique,  however,  permits  accurate  evaluation, 
as  follows. 


The  form  of  the  integrands  in  Eq.  (30)  is 

- fff  ( “x(->  ) 

C (p)'vJJ)  dQdrdp  rp  V e - 1 /•  J 


(krp/L  ) f j (r) 


C dedrdp  rp  J (krp/L  ) fj  (r)  f2  (p.p)  (76) 

*2  o 

where  fi(r)  - exp  | -r^-^-j  ~ “p2  O " L/F)2)  -2(r/oo)5^3}  (77) 

for  a general  focus,  and 

f2(P, p)  ' eik  p,°/L  H(F,p).  (78) 


Ue  note  that  S is  the  angle  between  p and  p. 

Each  of  the  integrands  contains  the  term 


rJ^f  krp/L  ) f j (r) 


Our  approach  is  to  expand  f^(r)  in  a Fourier-Bessel  series  and  then  utilize 


r J (or)  J (Sr)  dr  = 25(u-B) /(>+£) 

o o 


to  simplify  the  integration.  The  expansion  is  given  by 
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where 


The  utility  of  the  expansion  steins  from  the  fact 
ble  for  some  value  (A)  of  r that  is  readily  calcu 
the  other  integration  limits,  and  its  shape  is  su< 
sented  by  only  a few  (e.g.  6)  terms  in  the  series 
The  covariance  may  now  be  rewritten 

- f T ACy(LP  /k  A.p.e) 

CT(P)  * Vb„  \de  Le  X “ 


The  normalized  variance  is  simply 

, V*  T ACv<LCn/kA> 


1II-B.  Variance  of  Irradiance 


series.  The  strongly  saturated  case  can  be  handled  by  Eq.  (A5),  which  is 
readily  evaluated  using  standard  numerical  integration  techniques.  The  tran- 
sition region  between  the  unsaturated  and  saturated  conditions  can  in  prin- 
ciple be  handled  using  a general  expression  in  Kef.  17;  however,  this  func- 
tion is  highly  complicated. 

The  numerical  evaluations  of  the  variance  are  given  in  Section  IV, 
Figure  7. 


II1-C.  Cavariance  of  Irradiance 


The  unsaturated  two-source  log  amplitude  covariance  required  in  Eq. 
(83)  has  not  yet  proven  amenable  to  simple  series  expansion  representation. 

We  are  continuing  to  examine  this  problem;  in  the  meantime,  a two-fold  gaus- 
sian  quadrature  technique  is  being  used,  which  is  somewhat  time  consuming  or 
the  computer  but  which  is  resulting  in  numbers  that  agree  well  with  experi- 
ment (Figure  12,  Section  IV). 

The  saturated  case  is  further  complicated  by  the  fact  that  a two-source 
covariance  function  in  multiple  scattering  has  not  yet  been  derived  by  ana- 


Experimental  Results 


Introduction 


An  experimental  program  to  verify  the  theory  was  described  ir.  Kef.  1 
and  has  been  fully  underway  for  several  months.  The  measurements  are  made 
with  an  argon  laser  at  4880X  over  path  lengths  up  to  1 km.  The  facility  and 
instrumentation  are  as  described  previously,  with  certain  additions  noted 

below. 

In  the  course  of  these  experiments,  three  difficulties  have  been  un- 
covered. First,  it  has  been  found  that  the  argon  laser  has  sufficient 

% 

doppler-broadencd  spectral  linewidth  ('-10  GHc)  to  significant  ly  reduct  the 
variance  of  irradiance  observed  from  a diffuse  target  (speckles)  in  the 
absence  of  turbulence.  In  order  to  alleviate  this  difficulty,  the  laser  was 


u.  L.  Fried,  "Propagation  of  Spherical  Wave  in  a Turbulent  Meuium", 
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substantially  modified  to  include  the  insertion  of  a temperature  stabilized, 
intra-cavity  etalon  which  resulted  in  single-axial-mode  operation  and  elimi- 
nation of  this  effect. 

Second,  it  was  found  that  in  the  case  of  the  very  small  covariance 
scales  that  can  result  under  conditions  of  strong  turbulence  and/or  a large 
target  spot,  the  receiver  apertures  must  be  impractically  small  in  order  to 
resolve  the  covariance  curve  and  measure  the  proper  (unsmoothed)  variance. 
Similarly,  the  resolution  of  the  receiver  separation  must  be  very  high.  In 
order  to  make  such  small  apertures  and  highly  resolved  separations  practical, 
two  telescopes  were  fabricated  for  insertion  in  front  of  the  apertures.  This 
resulted  in  the  characteristics  shown  in  Table  III. 


TABLE  III.  RECEIVER  PARAMETERS  VS.  TELESCOPE  ATTACHMENT 


Effective  Receiver 
Magnification  Aperture  (mm) 

1 (no  tele-  2.0 

scope) 


Minimum 


Relative 


Separation  (mm)  Detector  Energy 


The  third  difficulty  is  more  fundamental.  It  is  of  course  desirable  to 
operate  in  all  six  parameter  realms  as  outlined  in  Table  I.  However,  as 
described  in  the  preceding  paragraph,  in  the  case  of  a large  target-spot,  the 
receiver  must  be  quite  small  and  will  in  principle  receive  a small  fraction 
of  the  transmitted  optical  power.  In  general,  if  the  target  spot-size  is  D, 
the  fraction  of  the  power  reflected  into  a Lambertian  hemisphere  which  is 
Intercepted  by  a receiver  of  size  L/kD  is  of  order 


L2/k2P2 

I* 


<v.  (target  spot-size  in  wavelengths) 


•«-  (kD)~2 


-2 


This  is  a small  number  and  independent  of  pathlength.  An  equivalent  viewpoint 

directly  relatable  to  Table  I is  that  a small  transmitter  of  size  (a^)  will 

result  in  speckles  of  the  same  size,  with  a fractional  receiver  energy  of 

order  (a  /L)  ; in  order  to  reach  parameter  realms  C,F,  and  E,  a will  be 
o . 2 ° 

constrained  to  be  smaller  than  / L/k  and  hence  (a  /L)  < 1 /L.  In  the  event 

o 

that  multiple  scattering  results  in  atmospheric  speckles  ((O  scalier  yet, 
the  problem  is  compounded  further. 

It  is  felt  that  the  current  experimental  program  is  sufficiently  com- 
prehensive to  verify  the  theory  presented  in  Section  II.  Certain  parameter 
realms  would  be  more  amenable  to  investigation  with  high-energy  pulsed 
lasers;  however,  the  realms  of  most  interest,  in  which  the  energy  is  concen- 
trated on  the  target  (focused)  to  whatever  degree  the  atmosphere  permits  in 
the  absence  of  an  adaptive  transmitter,  are  readily  handled  with  the  cw 
facility. 

In  addition  to  completing  the  experiments  described  here,  future  work 
will  involve  more  complex  targets  including  multiple  glints;  and  may  include 
wavelength  scaling  (e.g.  3.8y)  and  relative  translation  between  the  target 
and  receiver. 

IV-B.  Variance  Measurements 

Variance  measurements  prior  to  the  use  of  an  etalon  in  the  laser  are 
shown  in  Figure  6.  The  heuristic  curve  of  Eq.  (74)  is  also  shown,  withY; 
determined  from  Bmall-^  measurements  as  being  0.33.  The  value  of  Y 2 is 
taken  as  unity  for  simplicity;  it  will  begin  to  drop  from  that  value  at  the 
largest  values  of  the  abscissa,  where  p approaches  the  transmitter  size. 
Although  this  approach,  and  especially  the  use  of  a fixed  value  ofy;,  are 
somewhat  crude,  the  general  behavior  is  properly  predicted. 
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Values  obtained  over  a 500m  path  with  the  insertion  of  an  intracavity 
etalon  are  shown  in  Figure  7,  along  with  the  numerical  predictions  for  the 
unsaturated  and  saturated  cases  (Section  III).  In  the  transition  region, 
the  higher  of  the  two  curves  will  more  accurately  represent  the  true  com- 
posite curve,  within  some  limits  determined  by  their  respective  ranges  of 
applicability.  The  general  behavior  of  Figure  4 is  seen  to  be  verified. 

Note  that  the  simpler  theory  of  Reference  1 predicts  a variance  of  unity  in 
all  cases. 


IV-C.  Covariance  Measurements 


Typical  covariance  measurements  for  the  weak-turbulence  case  (param- 
eter realm  B-2,  Table  I)  are  given  for  two  different  pathlengths  in  Figures 
8 and  9.  The  theoretical  predictions  shown  are  from  the  joint-gaussian 
theory  of  Reference  1 and  in  particular  neglect  amplitude  perturbations. 

These  predictions  are  seen  to  be  accurate  except  for  the  tails,  as  may  be 
expected.  The  narrowing  of  the  curves  when  turbulence  increases  is  shown  in 
Figures  10  and  11  where  the  effect  of  decreasing  is  evident,  In  Figures 

12  and  13,  we  show  the  results  for  strong  turbulence  where  p is  less  than 

o 

and  comparable  to  /L/k  (realm  A-2  or  D,  Table  I).  In  this  case,  the 
serious  departure  from  the  simpler  theory  is  quite  evident.  The  substantial 
tail  on  the  measurements,  compared  to  the  joint  gaussian  theory,  is  explained 
by  the  general  theory  developed  in  the  present  report ; in  Figure  12  we  show 
several  points  obtained  with  the  numerical  techniques  described  in  Section 
III,  and  excellent  agreement  is  evident.  Note  that  the  (conventional) 
normalization  by  Cj2  obscures  the  further  deficiency  in  the  simpler  theory, 
i.e.  at  CjfO). 

In  Figure  14  we  show  the  measured  1/e^  points  of  the  covariance 
curves  as  a function  of  turbulence  strength  for  a pathlength  of  500m.  The 
theoretical  prediction  of  Reference  1 is  again  shown,  and  it  is  apparent 
that  the  results  in  strong  turbulence  are  determined  by  /l Tk  as  expected 
when  amplitude  perturbations  become  important. 

Finally,  in  Figure  15,  we  show  the  effects  on  the  tail  of  the  covar- 
iance of  arbitrarily  increasing  the  spot  size  (through  defocusing).  The 
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Figure  10.  Covariance  for  a focused  transmitter.  Moderate  Turbulence 
p • 2.6cm,  /L/k  c 6.2mm,  a ■ 1.35cm) 

O n ' 


'<■  iRurc  11.  Covariance  for  Focused  Transmitter,  Moderate  Turbulent' 


Theory-Joint  Gaussian  Assumption 


Experimental  Points 


Theoretical  Points,  Complete  Theory 


Covariance  for  Pot: used  Transmitter,  Strong  Turbulence 
P * 8.4mm,  /L/k  = 6.2mm,  a = 1.35cm) 


Figure  14.  Covariance  Scales  vs.  o 2 for  Focused  Transmitter  (I. 


figure  15.  Tail  of  Covariance  Funcli 


O 


Focused  Transmitter 
S 'v  1 cm 


S = 5 cm 
S = 10  cm 


vs . 


Increasing  Target  Spot  Size 


reduction  of  the  slope  and  level  of  the  covariance,  through  spot-size 
smoothing,  is  evident;  at  sufficiently  large  separations,  larger  spots  may 
be  expected  to  give  greater  residual  covariances  than  smaller  spots. 

IV-D.  Other  Properties  of  the  Irradiance 

As  discussed  in  Section  II,  wc  expect  the  probability  distribution  of 
irradiance  to  be  substantially  exponential  in  the  .ase  of  weak  turbulence. 

A typical  result  is  shown  in  Figure  16,  along  with  a plot  of  an  exponential 
distribution  having  the  measured  mean  value.  The  experimental  normalized 
variance  is  1.02,  while  unity  is  consistent  with  the  exponential  distribu- 
tion. 

In  stronger  turbulence,  we  expect  a combination  of  exponential  and 
log  normal  statistics.  A typical  experimental  result,  which  is  in  qualita- 
tive agreement  with  this  expectation,  is  shown  in  Figure  17.  Note  that  the 
values  at  small  irradiance  may  be  contaminated  by  system  noise.  The  experi- 
mental normalized  variance  for  this  case  is  1.2. 

The  existing  and  ongoing  data  will  also  be  used  for  the  computation 
o£  temporal  power  spectra  and  time-lagged  covariance  functions. 
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APPENDIX  A.  Generalized  Point-Source  Mutual  Coherence  Functions 

The  point-source  mutual  coherence  functions  are  crucial  to  the 
utilization  of  the  extended  lluygens-Fresnel  approach.  The  basic  work  on 

this  problem  was  done  by  Yura;  a complete  generalization  is  being  published 

19 

by  the  present  authors.  In  this  appendix,  we  summarize  the  second-oraer 
functions  of  interest  in  this  report. 

We  point  out  that,  following  the  completion  of  the  preceding  report 
(Ref.  1),  an  integration  error  was  discovered  in  Ref.  3,  and  the  expressions 
given  below  incorporate  appropriate  corrections. 

The  general  n -order  MCF  is 

H(pj  ,P2-  • -Pzni'Pl  1P2  • • -P?n) 

* |.«-<PJ,P1)+  p2>  + ...+  <!'(P2„-l.Pin-l)+  ^*(P2n.Pzn>;  > 

, (Al) 

■ where  1)  *=  x + !$• 

Under  certain  assumptions  discussed  in  Ref.  19,  the  primary  result  is 


2n 


H (p.p)  - exp<-  -r  E 
n 1 1 j=i+l 


2n-l 

l 

i=l 


K 20-2 


t (-l)1+3+1  D. 


+2  Z 


!■  C 


'ij 


j = l i = l 


<A2) 

where  M is  the  number  of  j satisfying  2j-l-i=cn  (for  i even)  and  2j+i=2n-l 
(for  i odd).  The  wave  structure  function  J,  and  log  amplitude  covariance 

■*  Q 

function  comprise  a generalization  of  those  normally  seen  in  the  litera- 
ture. They  refer  to  complex  phase  or  log  amplitude  quantities  as  a function 
of  (Pj-P^)  *n  t*le  observation  plane,  for  fields  emanating  from  and  in 
the  "source"  (target)  plane. 

The  phase  structure  functions  are  given  by 
1 


1/P  5/3  16/3 
0 


J <)t!(pi-pj 


)L  + (1-t)  (c.  J-C  j)  i 5/3 


(A3) 


19.  M.  H.  Lee,  J.  F.  Holmes,  and  J.  R.  Kerr,  "Generalized  Spherical  Wave 
Mutual  Coherence  Function",  to  be  published  in  J.  Opt.  Soc . Am. 
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where  p “ (0.545  C 2Lk2)  is  the  turbulence-induced  coherence  scale  and 
o n 

C 2 is  the  refractive  index  structure  constant, 
n 

The  log  amplitude  covariance  function  requires  some  discussion.  In 
the  first-order  or  vealc-scatterlng  case,  the  well-known  result  is 

1 oo 

CxiJ  - 0. 132w2k2L  Cn2J  dt|  du  u'8/3  sin2 
o o 


X Jo  [ul(Pi"P.j)t  + d-OCp^Pj)!  ] 


(A4) 


In  the  strong-scattering  (saturation)  case,  the  required  two-point  (p^,p.) 
function  has  not  been  derived.  However,  in  many  cases  the  one-point  func- 
tion suffices.  A modification  of  the  expression  given  in  Ref.  17,  such  that 
the  saturated  point-source  variance  is  unity,  gives 

1 “ I I 

.-3/5  |prv  > (A5) 

p„  ’ 


* OO 

((pi»Pj)  * 0.17^  du ^ dq  e q | 


5.69  uq 


This  expression  is  the  asymptotic  form  for  arbitrarily  strong  scatterings; 
a more  general  expression  is  given  in  Ref.  17.  The  latter  result  should  in 
principle  be  utilized  in  many  cases,  i.e.,  involving  multiple  scattering 
conditions.  However,  since  it  is  a relatively  recent  development,  it  has 
not  been  widely  used.  It  has  been  incorporated  in  the  present  report,  where 
appropriate.  Including  numerical  evaluations. 

We  now  consider  special  cases  of  interest: 

(a)  n - 1 

D 

H(Pi.P2;Pi.Pz)  “ e ’’12  (A6) 

0>)  n “ 2 

H^.Pa.Pa.P^Fj.pz.p^p,,)  - e"^(D12'Dn+Dl^+D2  3-D2^30+2C  +2C 

*13  *24 

t (A7 ) 

This  ia  a generalization  (double-argume.ics)  of  the  earlier  results  of 

Ref.  20. 

20.  D.  L.  Fried,  "Atmospheric  Modulation  Noise  in  an  Optical  Heterodyne 
Receiver",  IEEE  Trans,  on  Quantum  Electronics,  QE-3,  213,  June  1967. 
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where  v(t)  is  the  transverse  wind  velocity  along  the  path,  and  t is  the 
desired  time-difference  associated  with  (p^Pj)*  We  again  let  P;=P2«P3=P<- 
and  t.,  = tj3  = r.  The  special  cases  of  interest  are 
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BASE  UNITS: 

Quantity 


Unit 


SI  Symbol 


Formula 


length  metre 

mass  kilogram 

time  second 

electric  current  ampere 

thermodynamic  temperature  kelvin 

amount  of  substance  mole 

luminous  intensity  candela 

SUPPLEMENTARY  UNITS: 


A 

K 

mol 

cd 


plane  angle 

radian 

rad 

... 

solid  angle 

steradian 

sr 

... 

DERIVED  UNITS: 

Acceleration 

metre  per  second  squared 

m's 

activity  (of  a radioactive  sourcej 

disintegration  per  second 

(disintegration's 

angular  acceleration 

radian  per  second  squared 

rad/s 

angular  velocity 

radian  per  second 

rad's 

area 

square  metre 

m 

density 

kilogram  per  cubic,  metre 

F 

kg/m 

electric  capacitance 

farad 

A-«A' 

electrical  conductance 

siemens 

S 

A /V 

electric  field  strength 

volt  per  metre 

H 

V/m 

electric  inductance 

henry 

V-s'A 

electric  potential  difference 

volt 

V 

W/A 

rJcct/ir  resistance 

ohm 

V'A 

elcctromo!i\e  lorcc 

volt 

V 

W/A 

energy 

ioule 

J 

N-m 

entropy 

joule  per  kelvin 

N 

I'K 

force 

newton 

kgms 

frequency 

hertz 

Hz 

(cycle)’. 

illuminance 

lux 

lx 

1mm 

luminance 

candela  per  square  metre 

edm 

luminous  flux 

lumen 

Im 

cd-sr 

magnetic  field  strength 

ampere  per  metre 

VVb 

A'm 

magnetic  flux 

weber 

V-s 

magnetic  flux  density 

tesla 

T 

Wbttn 

magnetomotive  force 

ampere 

A 

power 

watt 

W 

)<* 

N m 

pressure 

pascal 

Pa 

quantity  of  electricity 

coulomb 

C 

A-s 

quantity  of  heat 

joule 

1 

N-m 

radiant  intensity 

watt  per  steradian 

Wrsr 

specific  heat 

joule  per  kilogram-kelvin 

I'kg-K 

stress 

pascal 

Pa 

Km 

thermal  conductivity 

watt  per  metre- kelvin 

W/m-K 

velocity 

metre  per  second 

ms 

viscosity,  dynamic 

pascal  second 

Pa-s 

viscosity,  kinematic 

square  metre  per  second 

ms 

voltage 

volt 

V 

W/A 

volume 

cubic  metre 

m 

wavenumber 

reciprocal  metre 

(wave)  m 

work 

joule 

1 

N-m 

SI  PREFIXES: 


Multiplication  Factors 


Prefix  SI  Symbol 


1 000  000  DOO  000  10” 
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• To  b*»  avoided  w|n?i#  possible 


tera 

T 

Riga 

C 

M 

mega 

kilo 

k 

hecto* 

h 
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da 

dm.l  * 

d 

cent  1 * 

c 

milli 

m 

micro 

M 

nano 

n 

pir  o 

fmnlo 

alto 

F 

a 
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MISSION 

of 

Rome  Air  Development  Center 


RADC  plans  and  conducts  research,  exploratory  and  advanced 
development  programs  in  command,  control,  and  communications 
( C 3)  activities,  and  in  the  C3  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  communications,  electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability , maintainability  and 
compatibili ty . 
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